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ABSTRACT 



A method for determining the free-fleld open-circuit voltage sensitivity 
of a reversible underwater electroacoustic transducer from the difference 
In Its input electrical Impedance when loaded by water and air was 
Investigated theoretically and experimentally. An equation for the 
sensitivity was derived using reciprocal two-port network theory. The 
theory takes into account the diffraction due to the finite size of the 
transducer, Its finite mechanical impedance, and Its free-fleld radiation 
impedance in water. An experiment to test the predictions of the theory 
was performed using a 6" diameter, hollow, piezoelectric ceramic spherical 
transducer. The result of the experiment agreed within several dB with a 
standard comparison calibration over the frequency range for which both the 
electrical impedance and comparison calibration data are considered 
reliable. The calibration method described, which has been termed the 
Delta-Z Method, may be useful for in-si tu monitoring of transducer 
sensitivity in Installations which can be flooded and purged. 
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I. INTRODUCTION 



A. BACKGROUND 

For the purpose of the research described in this thesis, an acoustical 
tranducer is a device which transforms fluctuations of environmental 
pressure into electrical quantities or vice-versa. A reversible acoustical 
transducer Is one which can act as either a transmitter or a receiver. 

There exist many calibration techniques which allow one to determine a 
transducer's receiving sensitivity, that is, the electrical output for a given 
pressure field at the transducer face, and/or Its transmitting response, that 
is, the pressure field resulting from a given electrical Input. The focus of 
this study is on a method Introduced by the Royal Australian Navy (RAN) 
whereby the open circuit voltage sensitivity of a reversible underwater 
transducer can be obtained from the difference In its Input electrical 
impedance (AZ) when loaded by water and air. This method shall hereafter 
be referred to as the Delta-Z method. 

While there Is agreement that the sensitivity can be obtained from AZ, 
there is disagreement over the appropriate analysis of the Delta-Z method 
and consequently over the formula which applies [Ref. 1] . Preliminary 
investigations of the Delta-Z method applied to in-service transducers were 
made by Taylor [Ref. 1] and Westbrook [Ref. 2]; however the experimental 
results were Inconclusive. 

Because of Its experimental simplicity, the Delta-Z method has great 
potential for in-situ transducer calibration. For this reason it was chosen 
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for further study. It was hoped that, by using a simple geometry transducer, 
for which many of the needed parameters can be easily calculated, the 
controversy over the proper analysis of the method could be resolved and Its 
usefulness established. 

B. OBJECTIVE 

The objective of the research described in this thesis is to develop, 
test and evaluate an underwater transducer calibration procedure for a 
simple geometry reversible transducer based on the difference In Its Input 
electrical impedance when loaded by water and air. 

C. OUTLINE OF REMAINDER OF THESIS 

In Chapter II the various standard methods used for transducer 
calibration and the theory associated with the Delta-Z method are 
discussed. The experimental procedure and the raw data gathered from the 
spherical transducer are described In Chapter III, while the analysis of this 
data and Its associated results are part of Chapter IV. A discussion of the 
results follows In Chapter V. A summary and some conclusions and 
recommendations appear in Chapter VI. 
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II. THEORY 



A. STANDARD CALIBRATION METHODS 

The calibration methods described in this section are discussed In 
greater detail by Bobber [Ref. 3] and some portions have been extracted from 
his book. 

I. Comparison Calibration method 

In this procedure the electrical output voltage of the transducer 
under test is compared of to that of a standard transducer for which the 
sensitivity is known. It is the most expedient and widely used method to 
determine the absolute sensitivity of most transducers. The results 
obtained using this method will be used throughout this thesis as the 
standard to which the results of other methods will be compared. 

The open-circuit output voltage, e s , of the standard hydrophone is 
measured in a free-field position with Its acoustic axis pointed towards the 
projector. The standard hydrophone is replaced by the unknown hydrophone 
and the open-circuit output voltage, e x , Is measured. If the free-field 
voltage sensitivity of the standard Is M s , then the sensitivity of the 
unknown is found from: 



Mx ~ ^s e x • 
e s 



(2.1) 
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2. Conventional Reciprocity Calibration Method 

A conventional reciprocity calibration requires three transducers 
of which one serves only as a projector (P), one is a reciprocal transducer 
(T) and serves as both a projector and a hydrophone, and one serves only as a 
hydrophone (H). Any one of the three transducers can be the unknown or the 
one being calibrated; however, the calibration formula usually is derived for 
the free-fleld voltage sensitivity, 1% of the hydrophone. The measurements 
are made in the far-field under free-f ield conditions. 

The arrangements and measurements made are shown schematically 
in Figure 1 and the procedure is as follows: 

a. the input current, ip^ to the projector is measured while the 
hydrophone free-fleld output voltage, ep(_|, is also measured , 

b. the same process is repeated with the hydrophone replaced by the 
transducer and its output voltage, epj, is measured, 

c. the projector is then replaced by the transducer and both the 
input current, ij, and the output voltage of the hydrophone, ejn» 
are measured. 
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Flgure I. Diagram of the Measurements Needed for 
a Conventional Reciprocity Calibration. 



The free-field open-circuit voltage sensitivity, M^, of the 
hydrophone is calculated using the following equation: 




e PH e TH j 



1/2 



e PT lT 



( 2 . 2 ) 



where J is the spherical reciprocity parameter: 

J= 2 A.- 
pf 

In this case, d is the distance between the acoustic centers of the 
two transducers, p is the density of the medium and f Is the frequency of 
operation. 
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3. i&oztraostiuceL R eci pro ci ty 

Only a single transmitting current and receiving voltage 
measurement (1j and ejH ) Is required to perform a reciprocity calibration 
If two Identical transducers are employed. In this case equation 2.2 is 
replaced by 



M H = 




1/2 



(2.3) 



B. SELF-RECI PROCITY METHODS 

1. Surface Reflected Self-reciprocity Method 

A single transducer can be used In a two-transducer reciprocity 
calibration procedure if It Is placed near a perfectly reflecting surface; 
reflecting a transmitted acoustic pulse back to the transducer, it receives 
Its own transmission. This self-reciprocity arrangement Is described by 
Bobber [Ref. 3] and Is shown schematically In Figure 2 where the air-water 
interface is used as the perfectly reflecting surface. 
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IMAGE 




Figure 2. Self-reciprocity Arrangement 



The resulting free-field open-circuit voltage sensitivity of the 
transducer, is determined by the expression: 



where: 




" e TH 



>T 



1/2 



ejH 3 open-circuit voltage across terminals during 
receipt of the reflected signal, 
ij = current flowing into the transducer during 
transmission, 



(2.4) 
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J s » spherical reciprocity parameter: 



J s =_2Ld_ (2.5) 

pf 

d = twice the distance from the radiating face of the 
transducer to the reflecting surface, 
p - the fluid mass density, 
f = frequency of the transmitted signal. 

2. Sabin’s Delta-Z Method 

Sabin [Ref. 4] pointed out that the ratio e^/it was equal to the 
difference the input electrical impedance before and after receipt of the 
reflected pulse: 

4Z.ZT-2f.em_, (2.6) 

It 

where 

Zj = input electrical impedance after receipt of the reflected pulse, 

If * input electrical impedance before receipt of the reflected 
pulse, i.e. in a free-field. The expression for the open-circuit 
sensitivity from equation 2.4 then becomes 

ri H * <AZJ S ) I/2 . (2.7) 
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3. Roval Australian Navy Method 

The Royal Australian Navy (RAN) developed a calibration method, 
claimed to be an extension of Sabin's method, In which the depth of water 
for the reflected Impedance measurement is taken to be zero (the 
transducer Is in air) [Ref. 5]. The planar reciprocity parameter Is used to 
determine the sensitivity instead of the spherical reciprocity parameter. 

The electrical impedance of the transducer to be calibrated is 
measured both in water (Z we t) and in air (Z^y) where: 

AZ = Z we t ~ Z(jpy , (2.8) 



so that the open-circuit sensitivity now becomes: 



M h = ( aZ J p ) 1/2 (2.9) 

where: 

Jp = the planar reciprocity parameter, 

Jp = 2A , (2.10) 

pc 

A - area of the active face of the transducer, 
p * density of the medium, 
c = speed of sound propagation in the medium, and 
M h * free-fleld voltage sensitivity. 

The method described above will be hereafter referred to as the 
Australian or RAN method, while the method that will be developed in the 
following section will be referred to as the Delta-Z method. 
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C. RECI PROCAL TWO-PORT NETWORK ANALYSIS OF THE DELTA-Z METHOD 
The treatment of the Delta-Z method as an extension of Sabin's method 
Is not justified. For one thing, in conventional reciprocity methods each 
transducer Is assumed to be In the far field of the other. Applied to 
reflection methods, this means that the transducer must be In the far field 
of Its Image In the reflecting surface. This condition cannot be met for a 
finite size transducer as the distance to that surface approaches zero. 
Furthermore, as will be apparent from the theory developed In this section 
and In the experimental results presented In Chapter V, the appearance of 
the planar reciprocity parameter in equation 2.9 Is Incorrect 
(notwithstanding that equation 2.9 does not Include the effects of 
diffraction and finite mechanical impedance). 

The proper analysis of the Delta-Z method Is In terms of a reciprocal 
two-port network. Such an analysis Is given In this section and leads to a 
different equation for the open-circuit sensitivity than equation 2.9. 

A reversible electroacoustic transducer may be represented as a 
reciprocal two-port network, as schown In figure 3: 
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Electrical 

Port 



Transducer 



Mechanical 

Port 




Figure 3. Two-Port Representation of an 
Electro Acoustic Transducer 

The appropriate two-port equations are 

e = Z E 1 ♦ T U, (2. 1 1 ) 

p = T 1 + Z m U, (2.12) 

here: 

e = the voltage across the electrical terminals, 

1 = the current passing through the electrical terminals, 

U ■ the volume velocity directed into the face of the transducer, 
p = the average acoustic pressure at the face of the transducer, 

\ * the electrical impedance with blocked mechanical terminals, 

Z m = the mechanical Impedance with blocked (open-circuit) electrical 
term Inals, and 

T= the transduction coefficient for a reciprocal transducer. 
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Now consider such a transducer radiating into air as schematized In 
Figure 4. 



Electrical Transducer 

Port 



Mechanical 
Port 



o 

+ 



U ♦— 



o 



e 



P=0 



o 



O 




Figure 4. Transducer Radiating In Air 



The radiation impedance of the transducer in air is considered to be 
negligible compared to Z m and hence Is represented in Figure 4 by a short 
circuit, 



p = 0 . (2.13) 

Solving for U In equation 2. 1 2 with p=0 and substituting In 2. 1 1 : 

U = zX±, 
z m 

e = Z E 1 + T (_£LL) . 

z m 
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Thus the electrical Input Impedance In air becomes: 



Zair^E'H • (2.14) 

z m 

Now consider the transducer transmitting Into a free-f leld in water, as 
represented by Figure 5. 



Electrical Transducer Mechanical Acoustic 

Port Port Field 




Figure 5. Transducer Radiating In Water 



Here Z ra( j ff is the free-f leld radiation Impedance of the transducer in 
water. Hence: 

p = -Z r adffU . (2.15) 



Substituting this Into equation 2.12 gives 



^rad ff U = T i + Z m U 
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or 



U- ill 



( z m + z rad ff) 



Substituting this into equation 2.1 1 yields 



e = Z E 1 + T 



_£Li 



z m + z rad ff 



The input Impedance therefore becomes: 



z water = -£- = Z E " 



z m + z rad ff 



(2.16) 



and the impedance difference between water and air is given by: 



AZ= z water " z air = Z E “ 1 2 



Z F + T 2 
z m + z rad ff z m 



which simplifies to 



1 2 z rad f f (2.17) 

^ z m + z rad ff ) z m 

Thus the electromechanical coupling coefficient T is proportional to the 
square root of the difference in the Input electrical impedance of the 
transducer when loaded by water and air. 
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Note that if Z ra( j ff and Z m are known, then T may be found from AZ 
using equation 2.17. If Z^ Is also known, then all the two-port parameters 

would be known and In principle one may analyze any proposed transducer 
application. In many cases of Interest, however, Z m and Z^ need not be 

known to obtain the sensitivity of the transducer, M 0 , provided that lAZI/IZI 

is significant. As will be shown below, this simplification is obtained when 
the radiation Impedance Is small compared to the mechanical Impedance, 
which is very often the case for piezoelectric ceramic transducers except 
near frequencies of mechanical resonance. This Is the great utility of the 
Delta-Z method. 

To relate the free field open circuit voltage sensitivity, M 0 , to AZ, 
imagine an acoustic wave Impinging on the transducer. The network 
representing this situation Is depicted In Figure 6. 



Electrical 

Port 



Transducer 



Mechanical Acoustic 

Port Fteld 




Figure 6. Transducer in an Acoustic Field 
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Where: 

Pb = the blocked pressure, that Is, the average pressure 
over the face of the transducer with zero 
volume velocity, 

pff = the free-field pressure, that is, the pressure which 
would exist at the location of the transducer if it 
were absent, 

D = the diffraction constant, and 
Z ra d ff = the free-field radiation Impedance of the 
transducer in water. 

By definition, the blocked pressure and the free-field pressure are 
related by the diffraction constant, D, which depends only on the geometry 
of the transducer and on the acoustic wavelength: 



Pb = Dpff. (2.18) 

The free-field open circuit (1=0) voltage sensitivity Is derived from 
Figure 6 as follows. 

The free-field open-circuit voltage sensitivity may be cast as 

M oc f f “ M o = ^oc_ - foe P Pb • (2.19) 

Pff P Pb Pff 

When i=0 is substituted in equations 2. 1 1 and 2. 1 2, they become: 
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eoc-™ 



giving: 



Also, from Figure 6, 



Therefore: 



giving: 



P = UZ 



m 



e oc =_L_ 
P z m 



p = Pb - Z rac jffU 



1 ♦ Z 



rad ff 



-m 



= Pb = Zm ♦ Z ra( j ff 
Zm 



Z 



m 



JL® . 

Pb Z m ♦ Z rad ff 



(2.20) 



( 2 . 21 ) 



( 2 . 22 ) 



Substituting equations 2.18, 2.20 and 2.22 Into 2.19, the open circuit free- 
field sensitivity becomes: 



T 



-m 



-m 



z m + z rad ff 



(2.23) 



The transduction coefficient T is found from equation 2.17: 



T = 



^Z Z m (Z m + Z rad ff) 



z rad ff 



1/2 



Equation 2.23 then becomes: 



”o = 



AZ z m (Z m ♦ Z ra( j ff) 



z rad ff 



1/2 D 



-m 



Zm (Z m + Z rad ff) 
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which simplifies to 



M 0 = D 



AZZ 






z rad ff + z rad ff) 



1/2 



(2.24) 



which is the desired result. 

If Z ra( j ff « Z m , as is usually the case for a piezoelectric transducer 

except near frequencies of mechanical resonance, equation 2.24 becomes 
much simpler: 



M 0 = D 



_AZ_ 



z rad ff 



1/2 



(2.25) 



Note that in this case the mechanical impedance need not be known to 
determine M 0 . This feature makes the Delta-Z method potentially very 
useful for in-si tu calibration of in-service transducers. Also note the 
similarity and difference between this and the Australian result for the 
sensitivity which was given in section II.B.3 and is reproduced here: 



M H = (AZ J p ) 1/2 



It is important to note that Z ra <j ff in equation 2.25 is not a reciprocity 
parameter and should not be construed as the inverse of the planar 
reciprocity parameter, J p , even if equation 2.9 took diffraction into 
account, which it does not. These differences will be apparent in Chapter V 
where the experimental results of both methods will be compared to the 
standard comparison calibration of a spherical transducer performed at 
TRANSDEC, San Diego. 
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a similar two-port analysis may be performed which relates the short- 
circuit free-field current sensitivity of a reversible transducer to the 
difference in the input electrical admittance when the transducer is loaded 
by air and water. The result is 



M 




AY Y rad ff 
( ,+Y m z rad ff)_ 



1/2 



(2.26) 



where: 

M sc = the short-circuit free-field sensitivity, 

AY = the change of admittance between the water and air, 

Y rad ff = the radiation admittance of the transducer, 

Y m - the mechanical admittance with shorted electrical 
terminals, and 

Ye = the electrical admittance with shorted mechanical 
terminals. 

This formula may be more appropriate for a low mechanical impedance 
device such as a moving coil transducer. (Note that a moving coil transducer 
Is anti -reciprocal and that we must replace T 2 in the above by -T 2 . This has 
no effect on the magnitude of the final result.) . 

In the case where 

Y m z rad ff < < ^ 



then equation 2.26 simplifies to 



M sc = D (AY Y rad ff ) 1/2 . (2.27) 
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The open-circuit voltage sensitivity can be determined using 
admittance parameters instead of impedance parameters with the following 
result: 



M 0 = zJH 

y E + (y E Y m~ j2)Z rad f f 



(2.28) 



where: 



T = the transduction coefficient for a 
reciprocal transducer, and 

T * t AY (V m *Y rad ff )] "2 . 
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III. DESCRIPTION OF EXPERIMENT AND RAW DATA 



A. DESCRIPTION OF THE TRANSDUCER 

The transducer used for this experiment was a 6-1nch outer diameter, 
hollow, piezoelectric sphere. It was borrowed from TRANSDEC and was 
built at NOSC, San Diego. It consists of two radially poled piezoelectric 
ceramic heml-spheres of unknown (probably about 3/8") thickness glued 
together. It Is potted in polyurethane of approximately 1/4 Inch thickness 
and has a 12 inch long hard rubber handle to house the connection to the 100 
foot long electrical cable (Cable * 71 102, DSS-3 (RF), MIL-C-915A), as 
depicted in Figure 7. Its low frequency capacitance is 210 nF (@1000 Hz). 

B. EXPERIMENTAL SET-UP AND PROCEDURES 

Measurements of Z wa t er and Z a j r were carried out at TRANSDEC, NOSC, 
San Diego. The main feature of this facility is a large water tank which 
simulates free-field conditions for a transducer placed at a certain position 
(6m depth at the center of the pool) [Ref. 6]. Its unique construction focuses 
the sound energy on the outside perimeter of the tank where it is reflected 
Into a sound trap, eliminating the energy reflected back to the transducer. 

Measurements were conducted using a continuous wave (CW) signal 
generated from the HP-4192 Impedance Analyzer which was controlled by a 
HP-9826 computer and HP Basic 4.0 software. The equipment set-up of 
Figure 7 was used to determine the resistance (R) and reactance (X) 
components of the impedance (Z) of the spherical transducer in both water 
and air. 
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Figure 7 : Equipment Configuration Used to Determine 
Sensitivity by the Delta-Z Method . 
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The procedure was as follows: 

The sphere was placed 6 m below the surface of the water for half an 
hour to allow its temperature to reach equilibrium with the surrounding 
water. The R wa ter and x water components were measured and stored on a 
floppy disk. The data covered the range of 1.0 kHz to 15.0 kHz in 100 equal 
increments. The sphere was then pulled out of the water, and R a | r and X a j r 
were immediately measured at the same frequencies as in water and stored 
on the same file as the water data. 

C. RAW DATA 

Figures 8 to 1 1 show the raw impedance data obtained for the sphere 
following the above procedure. 
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Figure 8. R W ater as a Function of Frequency. 
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Figure 9 . X wa t er as a Function of Frequency. 
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R air vs Frequency 




Frequency (kHz) 



Figure 10. R a j r as a Function of Frequency. 



X air vs Frequency 




Figure 1 1. X a j r as a Function of Frequency. 
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IV. AN ALY SIS ANP .BESULIS 



A. EXTRACTION OF AZ FROM RAW DATA 

The electrical Impedance measured at the terminals of the transducer 
Is represented by the complex quantity: 

Z = R ♦ j X 

where : 

Z = impedance, 

R = resistance, and 
X = reactance. 

For the measurements taken in water the expression Is written: 

^water = ^water * J x water 

and slmilarlly for the air measurements: 



^alr = R a1r + J x alr. 



The magnitude of AZ for Insertion In the sensitivity calculations is given 
by: 

I aZI = [(Rwater - Rair ) 2 + (Xwater -^alr^ 2 ^ (4.1) 

and the angle associated with It is: 
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IDelta-ZI/Z water 



(42) 



z AZ = arc tan 



*water ~ x alr 
^water ’ ^alr 



Only the magnitude of AZ was used in this experiment 

The values of resistance and reactance used in equation 4.1 were 
obtained as described in section 1 1 1.6. Figure 12 shows a plot of 
lAZl/IZwaterl vs Frequency. 



IDelta-Zt/Z water vs Frequency 




Figure 12. Plot of |AZI/|Z wa t er l vs Frequency 
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IDelta-ZI/Z water 



B. VARIABILITY OF AZ MEASUREMENTS 

To ascertain the reliability of the values of AZ, the experiment was 
repeated three different times. Figure 13 depicts the results obtained for 

lAZI/IZwaterl ^ rom tflese three independent experiments. 



LDelta-Z.l/Z water vs Frequency 




Frequency (kHz) 

Figure 13. Plot of lAZ|/|Z wa t :er l for Three Independent 
Measurements 

The results obtained indicate that random experimental error is 
probably less than 1% for frequencies below resonance. A detailed error 
analysis is included in Appendix A. Determination of the remaining 
parameters needed to calculate Mo (D, Z rac j ff, Z m ) is described below. 
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C. DETERMINATION OF D, Z rad ff and Z m 
1. Diffraction constant (D) 



RecallTrom Figure 6 and the discussion which followed that the 
diffraction constant Is defined to be the ratio of the blocked pressure to the 
free-field pressure. It was derived by Henriquez for a sphere [Ref. 7] and is 
given by : 

D = 1 (4.3) 

-\fT+ (kr) 2 

where: 

D = magnitude of the diffraction constant, 
k = wave number = «/c, 
a = radius of the sphere. 

The diffraction constant of a transducer with a more complex 
geometry is not necessarily a simple expression and may have to be found 
experimentally. Figure 14 depicts the diffraction constant as a function of 
frequency for several Interesting shapes. 
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Figure 14. Diffraction Constant as a Function of 
Frequency [Ref. 7] 

2. Radiation Impedance (Z^h tf) 

The radiation impedance of a sphere in radial motion is given by 

Beranek [Ref. 8]: 



where: 



Z rac j ff = Radiation impedance in free-field 
= Jcop 0 ca 
S (c + j«a( 1 -jcS)) 

O) 



(o = radial frequency, 
p 0 = density of the medium, 



(4.4) 
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c - speed of sound propagation in the medium, 

S = surface area of the sphere, 
a = radius of the spherical transducer, and 
- 6 - absorbtion constant. 

Since c8/<ois in the order of ixlO" 3 for all frequencies of interest in this 
study, it can be neglected without significantly affecting the results [Ref. 

8, pp 57,58]. Therefore it is assumed: 

Z ra d ff ~ j PqC (ka) ( 1 ~jka) . (4.5) 

S [ l+(ka) 2 ] 

The real and imaginary parts of equation 4.5 are : 

Rrad “ ppC(ka) 2 = (ka) Xp ac j (4.6) 

S(Wka) 2 ) 

X r ad = Po c Xa (4.7) 

S (l+(kr) 2 

3. Blocked Mechanical Impedance (Z m ) 

The blocked mechanical impedance of the sphere was extracted 
from its input electrical impedance in air and water as described below. 

The success of the technique employed depends upon the assumption that the 
radiation resistance in water is significant compared to the blocked 
mechanical resistance. 

The input electrical impedance of a simple reversible piezoelectric 
ceramic transducer in air may be modeled by the approximate equivalent 
circuit shown in Figure 15 [Ref. 9]. 
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Figure 15. Approximate Equivalent Circuit of a Simple 
Piezoelectric Transducer In Air. 

where: 

e 0 - open circuit voltage, 

C 0 = blocked electrical capacitance, 

Cn = blocked mechanical capacitance, 

I_n = blocked mechanical Inductance, and 



R^ = blocked mechanical resistance. 

The notation follows Wilson [Ref. 9], i.e. capital letter subscripts refer to 
equivalent electrical quantities, lower case subscripts refer to mechanical 
quantities. Only one mechanical resonance has been modeled In Figure 15; in 
principle any number of additional resonances can be modeled by adding one 
more parallel R-L-C branch for each. 

In making the above electrical analogy It Is assumed the vibrating 
elements are mechanically equivalent to a single degree-of-freedom mass- 
sprlng-dashpot system and that the radiation Impedance In air can be 
neglected. It Is also assumed that the equivalent circuit parameters are not 
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frequency dependent. The latter assumption Is not crucial to the success of 
the Delta-Z method, however, since the value of the mechanical Impedance 
significantly affects the calculated sensitivity only near frequencies of 
mechanical resonance. 

When the transducer Is surrounded by water the radiation Impedance 
cannot be neglected and Is represented In Figure 16 by and Iraq. 



Figure 16. Approximate Equivalent Circuit of a Simple 
Piezoelectric Transducer In Water. 

The procedure used to extract the mechanical Impedance of the 
sphere, or, more precisely, the mechanical resistance R m , Inertance L m , and 

compliance C m , Is as follows. 

At the frequency of mechanical resonance in air (« 0 ), the effects of 
and cancel and leave R^ in parallel with C 0 (see Figure 15), giving the 
following expression for the Input electrical impedance: 




Co 



3 



l l m + l rad 



6 



40 



(48) 



Zalr = R M ( 1 • 

I + (%R^C 0 ) 2 

where : 

«o = the natural mechanical resonance frequency. 

Separating Z 3 j r into its real and imaginary parts we obtain: 

^air = r m (4.9) 

I + (<o 0 Rfv 1 C 0 ) 2 

and X alr = -R h 2 fafrC 0 . (4 10) 

i + (%RmC 0 ) 2 

Rearranging equations 49 and 410 gives: 

X air s “tofyl c o . (411) 

R air 

R air 2 +x air 2 * R M 2 (4.12) 

1 + (« 0 R M c o )2 

Values of R a | r and X 3 j r values are obtained from the input electrical 
impedance data with the transducer in air. R^ and C 0 are then easily found 
by inverting equations 41 1 and 412. The same procedure is used to find the 
values of (R^j ♦ Rrad> and c o w ^h the transducer in water (Figure 16). The 
values of R^, C 0 , and Rp^o at resonance found this way are: 

R n - 1 29 n 
C 0 = 209. nF 

r rad = ^ J o 7 ^ 
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Using the fact that the transformation of the mechanical 
resistances, R m and R ra( j, Into their equivalent electrical quantities, R^ and 
Rrad> Is identical, the following may be written: 



R m = R M R rad • 1 3) 

r RAD 

Hence R m may be found from this equation using the values of Rr A q and R^ 
extracted previously once R rac j is known. Since R rac j can be easily 
calculated for the sphere, this is a simple procedure. The determination of 
R ra( j will be one of the crucial factors on which the successful application 
of the Delta-Z method to in-situ transducer calibration will depend. (It is 
not necessary to know R rac j, however to monitor degradation in sensitivity.) 

It should be noted that the ratio R m /R ra d * s related to the 
electroacoustic efficiency of a piezoelectric transducer, ii.as defined by 
Wilson [Ref. 9,p. 39] and Hunt [Ref. 10], by: 

R M s R m s 1 ~ *1 (4.14) 

r RAD R r ad *1 

so that equation 4.13 may be written 

R m = 1 ~ H R rad 

n 

Having found R m from equation 4.13, L m and C m are found from the 
center frequency, w 0 , and quality factor, Q, of the mechanical resonance 
using 

L m = Q R m (4.15) 

«o 
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and 



1 



(416) 



c m ‘ 

“o QR m 

The resulting values of R m , L m and C m for the sphere are: 

R m = 5.93 x 1 0 6 N sec/m 5 
L m = 1.25 x 10 3 Nsec/m 5 
C m = 1.46 x 10 -13 m 5 /N sec, 

and the electrical values of the electrical equivalent values R^, and 
are R^ ■ 1 29 0 

t M = 27.2 mH 
C(^j = 6.7 nF , 

Figure 1 7 shows a plot of Z m calculated using these values to the 
radiation Impedance calculated for the sphere. 
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Frequency (kHz) 



Figure 1 7 . Z m and Z ra(3 17 as a Function of Frequency 
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V. DISCUSSION 



A. COMPARISON OF DELTA-Z METHOD TO STANDARD CALIBRATION METHOD 

The standard of comparison for all other calibration methods used to 
determine the sensitivity of the spherical transducer is the comparison 
calibration described in Section II.A.1. The comparison calibration of the 
spherical transducer was obtained at the Transducer Evaluation Center 
(TRANSDEC) at the Naval Ocean Systems Center (NOSC), San Diego, Ca, using 
a Navy standard type F-37 transducer. The comparison calibration method Is 
widely used by the Navy and most other organizations that require absolute 
calibration of hydrophones. 

Figure 18 shows the results obtained at TRANSDEC by the Delta-Z 
method using eq. 2.24 (curve labeled "Delta-Z Method") compared to those 
obtained from the standard comparison calibration method. The two 
methods give results which agree within several dB except at the extremes 
of the frequency range of the experiment. The results from the Delta-Z 
method should be disregarded for frequencies above approximately 12 kHz 
since the second and third mechanical resonances occur In this region and 
were not accounted for In the estimation of Z m . Also, the results from the 
comparison calibration are suspect at frequencies below several kHz, as the 
pool at TRANSDEC cannot be considered anecholc at these low frequencies. 

In the frequency range 3 to 12 kHz the largest discrepancy between the 
results of the two methods is approximately 3 dB while the results over 
most of the frequency range agree within about 1 dB. 
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SENSITIVITY (dB ref lV/uPa) 



SENSITIVITY CALIBRATION CURVE 




FREQUENCY (kHz) 

Figure 18. Sensitivity Obtained From Various Calibration Methods 

The importance of including the effects of finite mechanical impedance 
in the Delta-Z method in order to obtain accurate results near frequencies 
of mechanical resonance is readily apparent. The curve labeled "Without Z m " 
in Figure 18 was calculated using equation 2.25, which assumes 
Z m >> Z r ad ff- The error made at 1 1 kHz by not accounting for finite Z m 
amounts to more than 5 dB. 
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The close agreement between the results of the Delta-Z and comparison 
calibration methods below the frequencies of mechanical resonance was 
quite unexpected and is very encouraging for the prospect of applying this 
method to the monitoring of in-service transducers. Before the Delta-Z 
experiment was performed it was anticipated that the precision of the 
results would be best for frequencies around mechanical resonance, as Z m is 
minimum there and the effect of water loading would therefore be more 
observable, i.e. lAZI/IZI would be larger. It was for this reason that a 
relatively large (by laboratory standards: 6“ diameter) sphere was chosen 
for the experiment. Based on the experimental results there is good reason 
to believe that if the diffraction constant and the free-field radiation 
impedance can be determined within a few percent, the sensitivity of an in- 
service transducer can be determined within 1 dB if AZ can be determined 
within 2 dB. From the magnitude of AZ/Z observed in this experiment 
(typically a few percent, see Figure 12), the precision with which 
impedance needs to be measured is probably about 0.1%. 

B. COMPARISON OF THE ROYAL AUSTRALIAN NAVY METHOD TO STANDARD 

CALIBRATION METHOD 

As was described in section II.B.3, the equation resulting from the 
Royal Australian Navy analysis of the Delta-Z method is: 



M|_| = ( AZ J p ) ,/2 



where: 



M|_| = open-circuit sensitivity of the transducer, 
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AZ = the change of input electrical impedance between 
water and air, and 

Jp = the planar reciprocity parameter. 

The results obtained using this equation are labeled "Australian Method" 
in Figure 18. The sensitivity calculated using the RAN method is in great 
disagreement with.that calculated using the other two methods. 
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VI. SUMMARY. CONCLUSION AND RECOMMENDATIONS 



A. INTRODUCTION 

The objectives of the research described In this thesis were 1) to 
develop an underwater transducer calibration procedure based on the 
difference in input electrical impedance when loaded by water and air, and 
2) to experimentally test the procedure using a simple geometry (spherical) 
transducer. 

B. SUMMARY 

Using reciprocal two-port network theory, an equation for the free-field 
open-circuit voltage sensitivity of a reversible underwater acoustic 
transducer was derived in terms of the difference in input electrical 
impedance when loaded by water and air. The theory takes into account the 
diffraction due to the finite size of the transducer, its finite mechanical 
impedance, and Its free-field radiation Impedance in water. An experiment 
to test the predictions of the theory was performed using a 6" diameter, 
hollow, piezoelectric ceramic spherical transducer. The results of the 
experiment agreed within several dB with a standard comparison calibration 
over the frequency range for which both the Delta-Z analysis and the 
comparison calibration data are considered reliable. 
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C. CONCLUSIONS 



It has been established that the open-circuit free-field voltage 
sensitivity of a reversible underwater piezoelectric acoustic transducer can 
be obtained by the so-called Delta-Z method, that is, from the difference in 
its input electrical impedance when loaded by water and air. The appropriate 
formula, derived from reciprocal two-port network analysis, is given by 



Mo = D 



A Z Z 



m- 



1/2 



z rad ff ^m + z rad ff) 



The finite mechanical Impedance of such a transducer must be taken into 
account in order to obtain precise results near frequencies of mechanical 
resonance. Apart from such frequency regions the radiation impedance is 
generally negligible compared to the mechanical impedance and the 
sensitivity can be found from the approximate formula: 



Mo = D 



AZ 

z rad ff 



1/2 



In this case neither the blocked impedance nor the open-circuit mechanical 
impedance need to be known to determine the open-circuit sensitivity. 

The Delta-Z method should not be regarded as a high precision method -- 
the precision of the results rests on the precision with which the difference 
in two quantities which often are nearly equal can be measured. 
Nevertheless, because it is such a simple procedure, it may be useful for In- 
situ monitoring of In-service transducers in installations which can be 
flooded and purged. 
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D. RECOMMENDATIONS 

Recommendations for further study are: 

1. Investigate the Delta-Z calibration method using piezoelectric ceramic 
transducers of more complex geometry with the aim of applying this 
method to the in-situ monitoring of in-service transducers. 

2. Investigate the application of the Delta-Z method to a low mechanical 
impedance transducer such as moving coil transducer 
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APPENDIX 



ERROR ANALYSIS 



A 6ENERAL 

Three types of errors are considered In this analysis: 

1. Random error, which is an indication of the repeatability of a 
measurement. This was determined experimentally by repeating the AZ 
measurement three times as described in Section IV.B. 

2. Estimated calibration error, which is an approximation of the 
systematic error of each measurement. 

3. Approximation error, values for the different parameters used in the 
calculation of sensitivity had to be truncated and/or estimated to a 
certain precision. 

B. ERRORINAZ 

lAZl “ [( R water" R air^ 2 + ^water"^air^ 2 ^ 

where: 

AZ = change of impedance for the transducer from water into air, 
R water = Resistance of the transducer in water, 

R a j r = Resistance of the transducer in air, 

X W ater = Reactance of the transducer in water, 

X a | r * Reactance of the transducer in air. 
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Using standard statistical analysis, the standard deviation for the aZ 
parameter is: 




V^water 



\6X a ir J 



The above assumes that R wa ten R air> x water and x alr are statistically 
Independent. 

If we further assume that 

So „ « Sd = Sy t = Sv . = S, then the fractional deviation of 

K water K a1r A water A air 

AZ is given by 



5 aZ _ 



AZ 



2 S 2 (R wa ter " R a1r^ 2 + 2 S 2 (X wa t e r ' x alr^ 2 
( R water ' R air^ 2 + ^water " ^air^ 2 



ni/2 



which leads to 
s az= 2S 



AZ 



AZ 



From the HP-4 192 A handbook; s - 0.3% of reading, so that 



SaZ " 0.6% 



AZ 
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C. ERROR IN THE BLOCKED MECHANICAL IMPEDANCE, Z m 



z m ~ R m + j( wL m + 1/o)C m) “ R m + J x m • 
where a is the radial frequency and L m and C m are given in terms of R m by 

L m = Q R m and 1 /( -m = w o ^ R m • 

« 0 

where: 

« 0 = mechanical resonance frequency (rad/sec) 

Q = quality factor of mechanical resonance. 

Then 

SL m ~ Q 5 R|yi and Stl/Cp^) = Wq Q fiRppj 
w o 

Ignoring the errors in f.a^ and Q, compared to the error in R m , 
then, 

SZpp| = I + j( 2nfQ + c 0 qQ) SR^ 

« 0 2nf 

and 

Sz m ~ SR m 
z m R m 
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Assuming an estimated 10% error on the blocked mechanical resistance R m 
(from the estimated 10% error on the equivalent electrical resistance R^), 

S7 ~ 10 % 

4 m 

z m 



It should be noted that the Influence of errors in Z m on the error in M 0 
will only be significant near resonance, where Z m Is of the same order of 
magnitude as Z rad ff. 



D. ERROR IN THE RADIATION IMPEDANCE, Z rad ff 

For the spherical transducer used In this project: 

Zrad ff = Pn c ka (ka + j) 

5 [1 ♦ (ka) 2 ] 

IZrad ff 1 * ppCka - (p 0 c) [ I * (ka)’ 2 ] ' m 

S [l*(ka) 2 ] l/2 S 

taking the log of both sides and differentiating: 

SIZrad ff ^ = ^Po^ ~ + 1 

lz rad ff 1 PqC S (ka) 3 [ 1 + (ka) -2 ] 

= p 0 c - .fii + _SCkal 1 
p 0 c S (ka) [ 1 + (ka) 2 ] 
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but = -2 = -2 , ignoring the error in k 

S a ka 



and 1 

[ 1 + (ka) 2 ] 



so that 



_$S. +Mal 1 

S (ka) [ 1 ♦ (ka) 2 ] 



lies between anc * 2 -^£2l 





ka 


ka 




then 








5z rad ff 4 


s 2 
b p 0 c 


♦ 2 fiflca) 2 


1/2 


z rad ff 


Po c 


ka 






2 s (ka). 


ignoring the error in p 0 c 




ka 







The fractional error in a, therefore In ka, is estimated to be about 5%. 
hence 



5|Z rad ff' - >0% 

*^rad ff 1 

E. ERROR IN THE DIFFRACTION CONSTANT, D 

D = I 
Vl * (ka) 2 
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Sq. (ka) 2 S (ka) - S( ka ) = 5.0% 
D 1 ♦(ka) 2 ka 
F. ERROR IN THE SENSITIVITY, M 0 

From Chapter II, 



M 0 = D 



AZ Z 



m 



_ z rad ff ^ z rad ff + z m). 



1/2 



For simplicity, define []= & z m 

z rad ff ^ z rad ff + z m) 

and z ra( j ff= Z r . 

Differenciating equation 2.24, 

S1 0 = Mo-ffi * ID [)-' n x 
D 2 

[]MZ + AZ6Z m - AZZ m (2Z(~SZ r + Z m 5Z r ♦ Z^Z m ) 
_ AZ Z r 2 +Z m Z r (Z r 2 + Z m Z r ) 2 

which can be written as 



5M 0 =5D + _D 8AZ ♦ Z r (5Z m /Z m ) ♦ (2Z r >Z m )(6Z r /Z r ) 
M 0 D 2 AZ (Z r +Z m ) 

but from sections D and E above, 



< MkSl and 6Z r N< 2 8(ka) 
D ka Z r ka 

so that 



(2.24) 
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8f = 1 + D (2Z r + Zp^) 5(ka) + D SAZ + Zp SZ^ 

2(Z r + Z m ) IT 2 AZ (Z r + Z m ) Z m 

Near resonance take Z r = Z m and D < 1. 



Then 




^Mo = [(2 x 5%) 2 + (0.5 x 0.4%) 2+ (0.5 x 10%) 2 ] 1/2 
Mo 



s Mo = 11% 
Mo 



Away from resonance let Z r « Z m and let (5Z m /Z m ) ~ (5Z r /Z r ) and take 
D« 1 , so that 

SMp = SD + D 5AZ ♦ 8Zr 
M 0 D 2 AZ Z r 

= 2 6(ka) + 1 SAZ 

ka 2AZ 



so 
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1/2 




S M 0 < 

M 0 



M 0 



'2 S(ka)\ 2 / Saz\ 2 
ka/ \2AZ> 

[ (2 x 5%) 2 + (0.5 x 0.4%) 2 ] 



10 % . 



1/2 



Converting the result to decibels: 

Error In dB = 20 log 1 +( S^ 0 /^ 0 ) = 20 log 1.1 = 0.8 dB 

Thus the estimated calibration error is approximately */. 1 dB for all 
frequencies. 
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